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Abstract 

(N 

^ \ New data on the anomalous magnetic moment of the muon together with the 

h — >■ Xg^ decay rate are considered within the supergravity inspired constrained 
. minimal supersymmetric model. We perform a global statistical analysis of 

\C \ these data and show that the allowed region of parameter space is bounded from 



below by the Higgs limit, which depends on the trilinear coupling and from above 
by the anomalous magnetic moment a^. The newest b — > data deviate 1.7 
a from recent SM calculations and prefer a similar parameter region as the 2.6 a 



p i' deviation from a^. 

• 1 Introduction 



' Recently a new measurement of the anomalous magnetic moment of the muon became 
I available, which suggests a possible 2.6 standard deviation from the Standard Model (SM) 
expectation Aa^ = a^^^ — = (43 ± 16) ■ 10~^°. The theoretical prediction depends 
on the uncertainties in the vacuum polarization and the light-by-light scattering, see 
e.g. the discussion in 0. However, even with a conservative estimate of the theoretical 
errors, one has a positive difference Aa^ of the order of the weak contribution to the 
anomalous magnetic moment, which opens a window for "new physics" . The most popular 
explanation is given in the framework of SUSY theories [0-|]l2|, since the contribution of 
superpartners to the anomalous magnetic moment of the muon is of the order of the weak 
contribution and allows to explain the desired difference Aa^. It requires the Higgs mixing 
parameter to be positive[|] and the sparticles contributing to the chargino-sneutrino (x^ — 
z/^) and neutralino-smuon (x° — /i) loop diagrams to be relatively light 0. 

The positive sign of /iq is also preferred by the branching ratio of the b-quark decaying 
radiatively into an s-quark - b —>■ Xg'y - ||13| . Last year the observed value of 6 — > Xg'j was 



close to the SM expectation, so in this case the sparticles contributing to the chargino- 
squark (x^ — q) and charged Higgs-squark [H^ — q) loops have to be rather heavy in 
order not to contribute to 6 ^ Xgj . 
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Figure 1: The dependence of the third generation Yukawa couphngs at the GUT scale 
as function of tan (3 for /xq > and /io < 0, obtained by fitting them to the low energy 
masses of the top, bottom and tau mass. The results are for a common mass mo = mi/2 = 
500 GeV, but for different masses the curves look very similar, except that the 'triple' 
unification point for /xq < shifts between 42 and 48, if the common mass is shifted from 
200 to 1000 GeV. 



However, it was recently suggested that in the theoretical calculation one should use 
the running c-quark mass in the ratio mc/mh, which reduces this ratio from 0.29 to 0.22 
g. The SM value for b X^^ increases from (3.35 ±0.30) x 10"^ to (3.73 ±0.30) x 10"^ 
in this case. This value is 1.7 a above the most recent world average of (2.96±0.46) x 10~^, 
which is the average from CLEO ((2.85 ± 0.35,^* ±0.22,j^,) x 10"^) ig, ALEPH ((3.11 ± 
QMstat ± 0.72,,,) X 10-^) and BELLE ((3.36 ± 0.53,t„t ± 0.42,,,(±ojO)„,^,;) x 10"^) 
r?! . For the error of the world average we added all errors in quadrature. 

As will be shown, the small deviations from the SM for both and h — > require 
now very similar mass spectra for the sparticles. 

In the Constrained Minimal Supersymmetric Model (CMSSM) with supergravity me- 
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diated breaking terms all sparticles masses are related by the usually assumed GUT scale 
boundary conditions of a common mass ttlq for the squarks and sleptons and a common 
mass mi/2 for the gauginos. The region of overlap in the GUT scale parameter space, 
where both and b Xg'y are within errors consistent with the data, is most eas- 
ily determined by a global statistical analysis, in which the GUT scale parameters are 
constrained to the low energy data by a minimization. 

In this paper we present such an analysis within the CMSSM assuming common scalar 
and gaugino masses and radiatively induced electroweak symmetry breaking. We use the 
full NLO renormalization group equations to calculate the low energy values of the gauge 
and Yukawa couplings and the one-loop RGE equations for the sparticle masses with 
decoupling of the contribution to the running of the coupling constants at threshold. For 
the Higgs potential we use the full 1-loop contribution of all particles and sparticles. For 
details we refer to previous publications [|18|, |T9[. 

In principle, one can also require b—r Yukawa coupling unification, which has a solution 
at low and high values of the ratio of vacuum expectation values of the neutral components 
of the two Higgs doublets, denoted tan/5 = {H^)/{H^)^ From Fig. |l] one observes 



that if the third generation Yukawa couplings at the GUT scale are constrained by the low 
energy top, bottom and tau masses, they become equal for /x < at tan/3 ~ 40, while for 
/i > they never become equal, although the difference between the Yukawa couplings is 
less than a factor three. Since /i > is required by Aa^ > (see below), we do not insist 
on Yukawa coupling unification and consider tan f3 to be a free parameter, except for the 
fact that the present Higgs limit of 113.5 GeV from LEPpO] requires tan/? > 4.3 in the 



CMSSM[13 



We found that the allowed area of overlap between b — > Xg'j and can be increased 
considerably for positive values of the common trilinear coupling Aq at the GUT scale. 
For Aq > the present Higgs limit becomes more stringent than for the no-scale models 
with Aq = 0, as will be shown. 



2 and b Xs^f in the CMSSM 

The contribution to the anomalous magnetic moment of the muon from SUSY particles 
are similar to that of the weak interactions after replacing the vector bosons by charginos 
and neutralinos. The total contribution to can be approximated by 



la^^^n tan/5 1 In """^ ~ 140 ■ 10"^^ tan/3 , 

^ Svrsm Owm%jjgY \ n J \ nisusY J 

(1) 

where is the muon mass, msusY is an average mass of supersymmetric particles in 
the loop (essentially the chargino mass). In our calculations we use the complete one- 
loop SUSY contributions from 0] with zero phase factors and the additional logarithmic 
suppression factor as in eq.(|T|). The calculated value of is shown in Fig. |^ as function 
of tan P . Clearly, it is approximately proportional to tan (3 and its sign depends on the 
sign of /ioB Only positive values of /io are allowed for the positive deviation from the SM 

^Our sign conventions are as in Ref. pl|. 
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Figure 2: The dependence of versus tan/5 for various values of the SUSY break- 

ing parameters mo and mi/2. The horizontal band shows the discrepancy between the 
experimental data and the SM estimate. Good agreement with the data is only achieved 
at large tan (3 and for light sparticles. Clearly, the fit allows only the positive sign of /x. 



and in addition the sparticles have to be rather light. However, light sparticles contribute 
also substantially to the b Xgj decay rate. In the past this posed a conflict. However, if 
one uses in the b Xg'y calculations the running mass for the charm quark, as suggested 
recently by Gambino and Misiak, the SM prediction is increased by 11%. In this case the 
newest world average on 6 — Xg'j is 1.7a below the SM, as mentioned in the introduction. 
Such a deviation is most easily obtained for large tan P and not too heavy sparticles, as 
shown in Fig. In the upper part the scale uncertainty of the low energy scale fib is 
displayed by the width of the theoretical curves, while in the lower part the dependence 
on the trilinear coupling Aq is shown. The scale fib was varied between O.Smf, and 2mb. 
For tan j3 ^ 40 only positive values of the Higgsmixing parameter at the GUT scale fiQ are 
allowed in agreement with the preferred sign of fiQ by the anomalous magnetic moment. 
For intermediate sparticle masses and /iq > large values of Aq and small values of the 
low energy scale {fib ~ CSm^) bring the calculated values of 6 — X^'j closest to the data, 
as can be seen from the left hand side of Fig. |^. Note that for heavy sparticles (right 
hand side of Fig. ^) the effect of the trilinear coupling is small, because the stop mixing 
is small, if the left and right handed stops are much heavier than the top mass. 
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Fig. ^ shows the values of 6 — Xg'y and a^^^^ as function of mo and mi/2 for tan/? =35. 
For 6 — >■ the ratio rndf^) / ml"'''^ = 0.22 was used, while for the NLO QCD contributions 
the formulae from Ref. ||2^ were used. The calculated values have to be compared with 
the experimental values BR{b X,7) = (2.96 ± 0.46) x 10"^ [HI- HI and Aa^ = 



(43 ± 16) ■ 10~^° |T[], which shows once more that b Xg'y and a^^^^ prefer a relatively 
light supersymmetric spectrum. 

To find out the allowed regions in the parameter space of the CMSSM, we fitted both 
the h — s> and data simultaneously. The fit includes the following constraints: i) 
the unification of the gauge couplings, ii) radiative elctroweak symmetry breaking, iii) 
the masses of the third generation particles, iv) h and Aa^, v) experimental 

limits on the SUSY masses, vi) the lightest superparticle (LSP) has to be neutral to be 
a viable candidate for dark matter. We do not impose h — t unification, since it prefers 
yUo < 0, as shown in Fig. |l], while Aa^ requires /iq > 0, as shown in Fig. ||. Yukawa 
unification for /io > can only be obtained by relaxed unification of the gauge couplings 



and non-universality of the soft terms in the Higgs sector p3 

The contributions of 6 — >■ X<j7 and the anomalous magnetic moment in the global 
fit are shown in Fig. ^for Aq = Q and tan j3 = 35. As expected, the contribution from 
b — > X<j7 is smallest for heavy sparticles, if 6 — X^'y is calculated with rric/mh = 0.29, 
while the minimum is obtained for intermediate sparticles, if rric/mh = 0.22 is used. 
With the newly calculated b Xg'y values, one can see, that b —>■ Xgj and prefer a 
similar region of the 1710,1711/2 plane. Fig. ^ shows the combined contributions from 
b — > Xs'j and a^^"^^ in the mo, mi/2 plane, both in 3D and 2D, for Aq = (top) and Aq 
free (bottom). In the latter case the lower 2a contour from b Xg'y moves to the lower 
left corner, but for the preferred value Aq ^ 3mo, which is the maximum allowed value in 
the fit in order to avoid negative stop- or Higgs masses and colour breaking minima, the 
Higgs bound moves up considerably. The total allowed region is similar in both cases, as 
shown by the light shaded areas in the contour plots. The 2a contours from the individual 
contributions are in good agreement with previous calculations 0, |[, but in these paper 
a simple scan over the parameter space was performed without calculating the combined 
probability. In addition, Aq = was assumed. 

We repeated the fits for tan (3 = 20 and 50, as shown in Fig. ^ For smaller values of 
tan P the allowed region decreases, since a^^ becomes too small. At larger tan /5 values 
the region allowed by and b —>■ Xg^ increases towards heavier sparticles, as expected 
from Eq. |l|, but it is cut by the region where the charged stau lepton becomes the Lightest 
Supersymmetric Particle (LSP), which is assumed to be stable and should be neutral. A 
charged stable LSP would have been observed by its electromagnetic interactions after 
being produced in the beginning of the universe. Furthermore, it would not be a candidate 
for dark matter. The increase of the LSP-excluded area is due to the larger mixing term 
between the left- and right handed staus at larger tan (3 . 

We conclude that the measurement strongly restricts the allowed region of the 
parameter space in the CMSSM, since it excludes the fiQ < solution, which was the 
preferred one from b — t Yukawa unification. In addition, it prefers large tan f3 with 
relatively light sparticles, if the present deviation from the SM of 2. 60" persists. 

At large tan /3 a global fit including both b — > Xg'y and as well as the present Higgs 
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limit of 113.5 GeV leaves a quite large region in the CMSSM parameter space. Here 
we left the trilinear coupling to be a free parameter, which affects both the Higgs limit 
constraint and the b Xg'j constraint, but in opposite ways, so that the preferred region 
is similar for the no-scale models with Aq = and models which leave Aq free. 

The 95% lower limit on mi/2 is 300 GeV (see Figs. 1^+0), which implies that the 
lightest chargino (neutralino) is above 240(120) GeV. The 95% upper limit on mi/2 is 
determined by the lower limit on a^^^^ and therefor depends on tan/3 (see Fig. For 
tan/? =35(50) one finds mi/2 < 610(720) GeV, which implies that the lightest chargino 
is below 500(590) GeV and the lightest neutralino is below 260(310) GeV. 
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Figure 3: The upper picture shows the dependence of the b — > Xg^y rate on tan/5 for 
Aq = and uiq = 600 (1000) GeV, mi/2 = 400 (lOOO)GeV at the left (right). For each 
value of tan /3 a fit was made to bring the predicted b — > X^'j rate (curved bands) as close 
as possible to the data (horizontal bands). The width of the predicted values shows the 
renormalization scale uncertainty from a scale variation between 0.5m;, and 2mb. 
The bottom picture shows the same dependence but for a fixed renormalization scale of 
Im^. The width of the band is given by the variation of Aq between -Stuq and Stuq. 
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Figure 4: The values of 6 — X<j7 and aj^ in the mo,mi/2 plane for positive /i and 
tan/5 = 35 to be compared with experimental data h Xg'-y = (2.96 ± 0.46) ■ 10^^ and 
a^^"^ = (43 ± 16) • 10~^°. One can see that both b Xs'j and a^^^^ prefer relatively 
light sparticles. 




Figure 5: The individual contributions to from b Xs'y and in the mo, Tni/2 plane 
for tan/5 = 35, > and Aq = 0. On the left handside we show the old contribution 
from b Xs'~f , as calculated with rric/mh = 0.29, which has the lowest fo^' heavy 
super symmetric particles. In the middle the contribution from b Xg'y for rric/mb = 0.22 
is shown, which now has a minimum for intermediate masses. The contribution from 
is shown on the right handside, which clearly prefers light sparticles. 
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Figure 6: The upper part shows the contribution (left) and its projection (right) 
in the mo,mi/2 plane for Ao = and tan/3 = 35. The light shaded area is the region, 
where the combined is below 4. The regions outside this shaded region are excluded 
at 95% C.L.. The white lines correspond to the "two-sigma" contours, i.e. = 4 for 
that particular contribution. The lower row shows the same for the fit, where Aq was left 
free, in which case Aq ^ Stuq (its maximum allowed value in our fit) is preferred in the 
region where the stop mixing is important, i.e. regions where the left and right handed 
stops are not very heavy compared with the top mass. One observes that with Aq as a 
free parameter the Higgs limit becomes the most important lower bound on the SUSY 
sparticles, while for the no-scale models with Aq = (top) the b Xg'-f rate determines 
mainly the lower bound. 
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Figure 7: The total ^ind the allowed regions in the parameter space for fi > and 
tan/5 = 20 (top) and 50 (bottom), with Aq free, as in Fig. ^ (bottom). 
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